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^ ' Cosmological perturbations due to statistical thermal fluctuations in a single fluid characterized by 

an arbitrary equation of state are computed. Formulas to predict the scalar and tensor perturbation 
spectra and nongaussianity parameters at a given temperature are derived. These results are relevant 
' to any cosmological scenario where cosmic structures may have been seeded thermaUy instead of 

originating purely from quantum vacuum fluctuations. 



I. INTRODUCTION 



o 

u 

In the standard inflationary paradigm, the seeds for cosmic structure are generated as quantum fluctuations. During 
Q ■ inflation, the quantum fluctuations of the fields present are stretched by the cosmic expansion to macroscopic sizes 
, and become classical [l|. These small inhomogcneities are then amplified in the later evolution of the universe 
C/3 ■ by gravitational collapse and eventually form the galaxies and other structures we observe around us today. The 
, ^ , \ predictions of simplest inflationary models can be matched with observations that require a nearly scale invariant but 
slightly red-tilted spectrum with only upper limits having been set on gravitational waves and any deviations from 
fSj \ the simplest statistical properties in terms of nongaussianity or statistical anisotropy Q • 

^ . Thermal fluctuations introduce another possible origin for small inhomogeneities and anisotropies. Thermal fluc- 
tuations are different from fluid hydrodynamical fluctuations [U, Q. In general, fluid fluctuations can arise from two 
different sources: There can be fluctuations in energy density and the associated temperature driven, for instance, by 
quantum fluctuations; this is what is traditionally discussed in literature. However, even if one can define a unique 
temperature in a given volume, there are fluctuations in energy within the volume due to the statistical nature of 
thermal physics. These are random fluctuations in all finite temperature systems that arise already at the classical 
\ level, and this is what is commonly referred to as thermal fluctuations ^ . In the early universe the temperatures could 
I . be very high, and therefore these fluctuations could be significant. The reason why in typical inflationary scenarios we 
• • ' do not worry about these fluctuations is that once inflation begins any "pre-inflationary" thermal matter is expected 
. ^ ■ dilute away rapidly leaving us with an almost pure vacuum state. However, there are cosmological models where 

■ thermal fluctuations could be solely or to a significant amount responsible for the initial seeds of inhomogeneities. 
j_] , For instance, in cyclic inflationary scenarios jJQ, where particle/entropy production keeps up with the inflationary 

■ dilution, thermal fluctuations become relevant ^ , for some of the interesting results we have found the reader is re- 
ferred to our companion paper Q. In bouncing cosmologies, where the big bang singularity is replaced by a smooth 
evolution from a contracting to an expanding phase, different matter sources become important near the bounce, for 
a recent review see 0], again making the thermal fluctuations relevant. 

In the present study we shall not refer to particular models but rather strive for generality. The purpose is to 
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develop the general formalism to tackle cosmological perturbations due to thermal fluctuations. We are going to make 
the following assumptions: 

• The universe contains a thermal fluid in "significant" abundance, i.e. the interactions within the fluid are able to 
maintain thermal equilibrium, this requires both kinetic and chemical equilibrium, for discussions, see [lol. [lll|. 
Typically this means that the relevant scattering rates have to be larger than the Hubble expansion rate. This 
requirement is often referred to as the ability of the fluid to maintain "local" thermal equilibrium. 

• For the sub-Hubble modes, (i.e. physical wavelength smaller than the Hubble radius, or the appropriate cosmo- 
logical time-scale), the statistical thermal fluctuations dominate over the quantum vacuum fluctuations inherent 
in the fluid. We will provide a quantitative criteria for this to occur in appendix IbI 

• There is no significant isocurvature perturbations due to the possible presence of other fluids. 

• There are no anisotropic stresses in any of the fluids. 

• There is some cosmological mechanism in place for the modes to exit from the sub-Hubble to the super-Hubble 
phase, after which the fluctuations evolve according to the usual hydrodynamical equations coupled to gravity. 
For instance, this could be attained by inflation where the Hubble radius stays approximately constant while the 
modes grow quasi-exponentially, or during a regular contraction phase where the Hubble radius shrinks faster 
than the contraction of the modes, as in the ekpyrotic [T2| or matter bounce Q scenario. For some other ways 
of realizing such a mechanism see hj. As previous literature dealing with thermal fluctuations 

EMS, 

we also 

assume that the transition from sub- to super-Hubble phase is instantaneous. To study the precise nature of 
the transition would involve non-equilibrium thermodynamics in curved space-time which is clearly out of the 
scope of the present paper, but we do not expect the results to be affected beyond 0(1) factors. 

• At least, near the sub to super transition we can trust General Relativity and the usual laws of thermodynamics. 

In previous literature thermal statistical fluctuations have been considered in a variety of contexts: The earliest 
hint that these could be relevant for CMB was perhaps provided by Peebles [11] , and further developed in [l^ [iGl lTgl 
(see also [131). Since then thermal fluctuations have found applications in several models: string cosmology |l5l. l2ll| 
inflation inflation (in particular warm inflation 



glUll - for reviews see e.g. [1,13 - but see also j23-l27j). bouncing 



cosmologies [T3| Milne/holographic universe In the present paper we first generalize the calculation of the 

curvature perturbations to the case when the thermal matter could have an arbitrary equation of state in the 
process clarifying several conceptual issues related to gauge choices, and transfer of perturbations from sub- to super- 
Hubble phase. These generalizations can be particularly important and interesting for early universe cosmology where 
stringy thermodynamics jsil - lssj and/or phase transitions may be relevant, neither of which is described by a constant 
equation of state parameter which is what previous studies have been mostly confined to. Next, we provide general 
formulas to compute not only the scalar power spectrum, but also the spectrum of gravity waves and higher point 
correlation functions. These results when applied to phase transitions in cyclic inflationary models turn out to produce 
interesting signatures for Planck, and will be discussed in Q. 

The paper is organized as follows. In section |TT] we will derive the curvature perturbation in a universe filled with 
a thermal fiuid with an arbitrary equation of state. We also compute the spectrum of gravity waves expected in this 
set-up, or equivalently the tensor-to-scalar ratio. These derivations are somewhat technical, but only familiarity with 
standard cosmological perturbation theory is assumed^. In section iHll we will derive the nongaussianity parameters 
due to the thermal fluctuations, and in section lTVl we illustrate the application of our formulas for a radiation dominated 
contracting universe. Section |V] briefly concludes. Appendix |^ concerns a technical issue of going over from real to 
Fourier space that is needed to make contact with the usual cosmological perturbation analysis; in appendix [B] 
we compare the relative strengths between quantum/hydrodynamical and statistical thermal fiuctuations, and in 
appendix [C] we calculate thermal pressure fluctuations for completeness. 



^ This is going to be particularly relevant for applications to cyclic inflation models. 

* For a very pedagogical and transparent introduction, see theory.physics.helsinki.fi/'^genrel/CosPerShort.pdf 



II. THE CURVATURE PERTURBATION FROM THERMAL FLUCTUATIONS 



A. Curvature Perturbation &: the appropriate Gauge Choice 



We are going to consider a cosmological set-up where the dominant fluid component of the universe is thermal, 
i.e. there exists local thermal equilibrium, so that as long as the wavelengths of fluctuations are smaller than the 
cosmological time scale their power spectrum is determined by the thermal fluctuations in the thermal fluid. Once the 
modes become super-Hubble, thermal correlations over the relevant physical wavelengths can no longer be maintained, 
instead the fluctuations evolve according to the usual hydrodynamical differential equations coupling the metric and 
the matter fluctuations. Essentially, in this set-up the thermal fluctuations act as initial conditions to seed the 
supcr-Hubblc fluctuations. 

Now, the supcr-Hubblc modes are easy to track because they behave as zero modes and it is well known [l| that 
the curvature perturbation, Cfc, remains a constant even if the equation of state parameter is not ^. In fact, the 
above statement is true even if General Relativity is not valid [42|, but as long as we are only looking at adiabatic 
super-Hubble perturbations. This makes our analysis applicable to several bouncing/cyclic models which resort to 
modifying gravity to obtain a non-singular bounce (modulo the caveats about mode mixing mentioned in the previous 
footnote). However, if the universe docs contain more than one type of fluid/flclds, then isocurvature perturbations 
could be important, but we are not going to consider them in this paper. 

Our goal in this section will be to compute (^k arising from thermal fluctuations in the sub-Hubble phase where 
thermal correlations can exist. In particular we will evaluate this at the Hubble crossing which, according to our 
previous discussion, will provide us with the primordial spectrum for CMBR. In the next section, we are going to 
calculate the two and the three point correlation functions as well as the gravity wave spectrum. We do not make any 
assumptions about whether we have an expanding or a contracting universe (again modulo the previous comments 
about mode matching) or whether there is a single thermal fluid or several energy components. However, we shall 
make the crucial assumption that the fluctuations are dominated by the thermal fluctuations of a single fluid, dp, 
where p is the energy density of the thermal fluid. We parameterize the contribution of the thermal fluid to the energy 
budget by 

where a is the scale factor of the universe, H = a/a is the conformal Hubble rate, Mp is the reduced Planck mass, 
and here and in the following the overdot denotes the derivative with respect to conformal time. Thus, if the thermal 
fluid is the only component in the universe, f2 = 1. 

We would now like to relate the perturbations in the fluid to a gauge invariant degree of freedom describing the 
metric perturbations. The scalar perturbations in the metric can be parametcrised as [l| 

ds^ = a^{ri) [-(1 + 2(j))dTj^ + B^.drjdx' + ((1 - 2i;)6,j + E^,j) dx'-dx^] . (2) 

We need also to parameterize the perturbations in the matter content, which will be treated as a perfect fluid, so no 
anisotropic stresses are present. Then the energy momentum tensor can be written as 

T%^~p{l + d), T'',=p{l + w)v.,, T)^p{w + clS). (3) 

Since we assume a thermal fluid, all the background quantities are given solely by the temperature, 

p = p{T), p = p{T). (4) 

and so is w(T). and c^(T') ^. In the following the prime will always denote the derivative with respect to the 
temperature, and we often drop the explicit argument T. An overdot will refer to a derivative with respect to the 
conformal time rj. 



In a contracting universe one of the two modes of (Jj. is growing (the one which is decaying in an expanding universe) , while the second 
is constant l39l |40|| . We will here assume that the growing mode in the contracting phase couples only to the decaying mode in the 
expanding phase. Whether this is the case or not depends on the specific model of the transition between contraction and expansion. 
For a discussion of this issue see e.g. |41| . 

Unlike in the case of hydrodynamical fiuctuations, the pressure fluctuation is not given by the adiabatic value, 5p = [p' (T) / p' {T))5p, 
but is rather determined via stress-energy conservation, see appendix \C\ for the precise formula. We will see shortly, that to compute 
the spectra we do not need the explicit form of the sound speed squared. However, as the modes exit the horizon, the perturbations are 
expected to relax to their adiabatic value which guarantees the constancy of the curvature perturbation at large scales. 
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Now wc want to calculate the gauge invariant curvature perturbation, which in a general gauge reads as 

C = - Hiv - B) . (5) 

Since it is a gauge invariant quantity we can evaluate it in any gauge of our choice. Let us choose to work in the 
longitudinal gauge where E = B = 0. In this gauge the metric can be written in terms of the gauge invariant Bardeen 
potentials $ and 



ds^ = a^ii-i) [-(1 + 2^)dr]^ + (1 - 2^)dx^] , 
and the Oi component of the Einstein's equations determines v in terms of the Bardeen potentials: 



(6) 



2M2 



(1 + ■w)pv . 



Using this in (O we have that 



2M^n 

(1 + iy)a2p 



(7) 



(8) 



Please note that the above equation is written in a completely gauge invariant form and is therefore valid in any 
gauge. This is important for us because it will become necessary for us to switch to the comoving gauge on physical 
grounds. This is actually a subtle issue which, to our knowledge, has not been explained before. The point is that all 
thermodynamic calculations, such as those relevant when we will derive the energy fluctuations in a given volume, are 
typically carried out in Minkowski space-time. In order to generalize the analysis to the FLRW metric (or any other 
metric for that matter) one has go to a frame where the background fluid is at "rest" [43| . which is none other than 
the comoving gauge. In this gauge Sp becomes gauge- invariant, as it must because the Minkowski calculations cannot 
obviously depend on the gauge choice, and is related to the Bardeen potential via the relativistic Poisson equation: 



1 



= — 



2 XkAL 



(9) 



The superscript "C" refers the comoving gauge which we are going to subsequently drop as all the thermodynamic 
calculations implicitly assumes this same gauge choice for the perturbation in the matter density field. 

At this point it is useful to set $ = 4', since we assumed that the anisotropic stresses can be neglected. We can 
then compute hi terms of the density fluctuations from ^ and substitute it in ([5]) to obtain 



c 



1 



2 \kMr: 



1 



{l + w)p 



(3 + r) 



Sp, 



where time evolution of the density fluctuation is parameterized as 

_ dlogSp _ {Spy T 
dloga Sp H 

We remind the readers that the prime corresponds to the derivative with respect to the temperature. 
More succinctly. 



(10) 



(11) 



Spk , 



where we have defined a time/temperature dependent proportionality coefficient 



2(3 



3{l + w)n 



(12) 



(13) 



for later convenience. These quantities will depend on the temperature at the time of the "exit" of a given comoving 
mode. 
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B. Thermal Density Fluctuations 



We will now use the thermodynamics to quantify fluctuations in the fluid and then use the results of the previous 
section to relate them to the metric perturbation spectra. This is more of a review of what has been discussed in the 
previous literature [isl - fTTj . and our results agree to within 0(1) factors, until a crucial step highlighted in the end 
this subsection. 

One defines the average fluctuation in energy, /S.E via 

{AE), = {E)-{E)) _ _^^_-^_TC^, (14) 

^^P)^ = ^ = T^df^ (15) 

where Cl is the heat capacity of the thermal system for a given volume . Note we have also introduced a subscript 
L in [lS.E)\ to denote that we are considering fluctuations in a given volume. 

The next step is to go from real space to momentum space. This is a tricky procedure as it depends to some extent 
on the window function one chooses. In appendix \^ we consider in details this procedure using gaussian window 
function and obtain 



5pl = p ( V>L=a/fc with 7 = 2y27r3/4 « 6.7 . (16) 



Thus we have 



leading to 



5pl ^ ^T^p' (17) 



and eventually 



using the standard definitions of the power spectrum. The subscript k refers to the fact that all these quantities have 
to be evaluated at the Hubble crossing condition, Hk ^ k/a, which we have also used along with ([T]) to eliminate the 
Hubble factors. A few comments arc now in order. Firstly, as in all previous literature on the subject in deriving 
(|19p we have implicitly assumed an instantaneous transition from the thermally correlated sub-Hubble phase to the 
hydrodynamical super-Hubble phase. This is obviously not realistic but a careful investigation of such a transition 
is very challenging because it will involve non-equilibrium thermodynamics on curved space-time, clearly out of the 
scope of the present paper. 

Secondly, one can see that the factor A{T) basically tells us what is the difference between the spectra of the 
gravitational potential (or, via the Poisson equation, the density perturbation) and the spectra of the gauge invariant 
curvature perturbation. It is crucial to take this into account: one can imagine physical situations where A can 
even vanish or diverge. This is where our results differ significantly from previous literature which typically only 
computes $fc and are directly applicable only for constant equation of state parameters when making comparison 
with observations. 

C. The prefactor A{T) 

The last missing piece required to obtain the power spectrum is an expression for A{T). Explicitly, our definition is 

6(1 + iu)U 



^ Different schemes typically yield slightly different values for 7. 
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w can be obtained straight forwardly as a function of temperature from the partition function, or p{T). A useful 
thermodynamic relation, in this context is 



so that 



w 



P 



Tp' — p 



(21) 



(22) 



The computation of Q at the exit temperature depends on the specific model under consideration, and one cannot 
make any further simplifications at this point. Obviously, if the thermal fluid is the only energy component in the 
universe, then = 1. 

We arc finally left with the evaluation of r. From the expression we obtained for thermal energy density fiuctuations 
(dZl) we first find 



Now, recalling the continuity equation 



2p' + Tp"\ dhiT 



2p' J dln( 



(23) 



we see 



Thus we finally have 



c^lnT p 

-3(1 + ^0^. 



din a 



1 + 



{l + w)p{2p' + Tp") 
Tp'^ 



(24) 



(25) 



Note that the sign of this remains the same also in the contracting phase: then the temperature is getting lower with 
time, but also the Hubble rate is negative. We will illustrate the computation of the power-spectrum for the special 
case of radiation towards the end of the next section. 



D. Gravity Waves 



Another interesting probe of our early universe is the primordial gravitational wave which are stretched like the 
scalar perturbations. However, for a linearized Einstein's gravity there is no source term for the gravitational waves. In 
principle the initial conditions for the gravitational waves could be set purely classically |44| or from quantum vacuum 
condition [l| ® . Assuming that the initial conditions for the primordial gravitational waves are set by quantum vacuum, 
i.e. Bunch-Davis, the gravitational wave spectrum is given by 



1 

4^ 



H 



The tensor to scalar ratio is then given by 



Vq 7' 12V3^3j^iA2(T) MpT^p' 



(26) 



(27) 



In general, the temperature dependence of this and the scalar spectra too depend very nonlinearly on the properties 
of the thermal fluid, but these are straightforward to compute once we know p{T). 



In fact, any thermal matter can only act as sources of gravitational waves if it's partition function is non-extensive [ISl. l32l|. a scenario 
not considered here. 
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III. NONGAUSSIANITY: BI-AND TRI-SPECTRUM FOR THE CURVATURE PERTURBATIONS 



In the previous section we evaluated the CMB power spectrum as a two step process. In section Hi Bl we calculated 
the thermal density fluctuations from the partition function (or pressure) governing the thermodynamics of the fluid 
in the comoving gauge in which the fluid is "at rest" and therefore the Minkowski space-time calculations can be 
applied [1^. In section Hi Al we found how these thermal fluctuations are related to the curvature fluctuations, which 
then allowed us to obtain the two-point correlation function in CMB. We can apply the same prescription to obtain 
higher point correlation functions, we just have to compute the appropriate higher thermodynamic cumulants (see 
also [IJI for an earlier study of non-Gaussianities from thermal fluctuations).. 

The third and the fourth order centered cumulants are given by 



In Z 



54 \nZ 



= {E-^)-2,{E^){E)+2{EY = {l\E-^), 



{E*) - A{E^) (E) + Q{E^){EY - A{E) {Ef = {AE^^ 



(28) 
(29) 



From the above thermodynamics we infer that 



(2/7' + Tp") 
2T^ (3p' + 2,Tp" + p'") 



(30) 
(31) 



where one considers fluctuations in a box of size L. These formulas hold in the real space, but can be converted to 
momentum space using window functions as before ^ we have 



(32) 



Using the standard definitions of the spectrum and nongaussianity parameters, see 0, |4^, we are ready to write 
down the results using in ([5^ : 



Jnl = 



9NL 



5fc-i (a) 



54 



1 



2\2 



1 



5p{2p' + Tp") 
2AT{p'f 
25p2 [3(p' 



(^2)3 1727^2 (T) 



Tp") 



T^p" 



243T2 {p'Y 



_ G{T) 



1727^2 (T) 



(33) 
(34) 



In deriving the expressions for the nongaussianity parameters we have used (jl2p to relate fluctuations in energy density 
to the curvature fluctuations, and also ([T]) to eliminate the Hubble factors. 

Physically, the most important aspect about the thermal NG parameters is that they depend on how the pres- 
sure/density varies as a function of the temperature. Moreover, the higher the order of the correlation function, the 
higher the derivatives that becomes relevant. It's not hard to realize then that if the exit temperature is close to a 
thermal phase transition, we might be able to see it's effects in the enhancement of NG parameters. This is precisely 
what we observed in the context of cyclic inflation scenario 0. For the present, we will just illustrate the various 
computations of the cosmological observables by considering the ordinary relativistic fluid. In particular, we will see 
that for pure radiation, the nongaussianities and the tensor-to-scalar ratio are too small to be observed by the Planck 
experiment. 



® The relations to the Fourier space spectra are tricky, and the 7 factors for the different correlations functions could be different depending 
upon the window functions used, but the difference is only expected to provide 0{1) modulations. 
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IV. EXAMPLE: RADIATION 



For the purpose of illustration let us consider a radiation dominated contracting universe. Since in this case the 
Hubble radius contracts faster (as 1/t) than the physical wavelengths (as the latter is pushed out of the Hubble 

radius and the various perturbative spectra becomes imprinted at the time of the mode-exit. We can use the formulas 
of the previous section to compute the different cosmological observables. 

For relativistic radiation fluid we have that 

p(T)=.gT4andp(T) = |r4. (35) 

where g depends on the number of relativistic degrees of freedom. It follows immediately that 

w = i and r = -4 , ^ A = i . (36) 
3 4 

In general, these parameters need not be constant, but they happen to be in this simple case, or whenever pressure 
is a power-law in temperature. 

The amplitude of the primordial spectrum, according to our convention, is then given by 

where T corresponds to the temperature when the given mode exits becomes super-Hubble. Evidently, the spectrum 
is not scale invariant because the amplitude depends strongly on the temperature, and T cx 1/a giving rise to a 
very large blue tilt. We should point out that our claim that ((37)) is the primordial spectrum relevant for CMBR 
relies crucially on the fact that there is no mixing between the mode of C,k which is growing in the contracting phase 
with the dominant mode in the expanding phase, the constant mode, and in addition, on the assumption that there 
are no entropy modes which become important and which could change the spectrum of the curvature fluctuations 
on super-Hubble scales. If there is unsuppressed mixing between the growing mode in the contracting phase and 
the constant mode in the expanding phase (see [46j for examples where this is the case), then the amplitude of the 
resulting curvature fluctuations changes, but the spectrum remains as given by ([37]). The reason that there is no 
change in the shape of the spectrum (in contrast to the case of a matter-dominated phase of contraction where the 
index of the power spectrum changes by —2, see [1^, E^) comes from the fact that for a radiative equation of state the 
canonical fluctuation variable v [43, which is related to C via C, ^ a~^v has vanishing squeezing factor and hence 
remains conserved. Thus, there is no preferential growth of long wavelength modes compared to short wavelength 
modes which would come from the fact that long wavelength modes spend more time on super-Hubble scales in the 
contracting phase. For a discussion of this point see [i^ . 

For radiation we obtain the following numbers for the nongaussianity parameters: 

fNL = — «0.16, (38) 
247 

25 

ff^L = ^ « 0.003 . (39) 

Not surprisingly for radiation, which is free from any internal scales, both the /jvl a-nd the g^vL parameters turn out 
to be scale invariant. The above approximate values correspond to the 7 value for the gaussian window function (see 
the appendix [^ which is unfortunately beyond Planck's sensitivity. 

Let us compute the tensor-to-scalar ratio for radiation. We readily obtain 

For a given temperature, we can fix the unknown g by matching the amplitude of perturbations p7p with the observed 
one. This then fixes the tensor-to-scalar ratio (|40p . In other words, we can deduce the primordial temperature from 
observations: 



In the second equality we have used the gaussian window value for 7, and the lower bound is obtained from the present 
best-fit WMAP value for the amplitude Aq = 2.4 • 10~^ and the bound rt/s < 0.24, which both apply at the scale 
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k = 0.002/MpC. The minimal allowed temperature corresponds to a huge number of effective degrees of freedom, 
g ~ 10^^. For p ~ 1 of order unity, the observed amplitude requires T/Mp ~ 10"'', in which case the tensor-to-scalar 
ratio will be too low to be measured. 

To conclude, thermal fluctuations in usual radiation cannot account for the CMBR spectrum as the spectrum 
is heavily tilted, not surprisingly. We also found that such fluctuations cannot produce large nongaussinities or 
gravity wave signals. Similar conclusions hold for simple polytropic thermal fluids, but as we will see in richer 
thermodynamics may indeed yield detectable nongaussianities. 



We considered statistical thermal fluctuations as a possible source for cosmological large scale structures. We 
presented a robust derivation of scalar and tensor spectra in this context. We also explicitly provided the formulas 
for the bi- and trispectrum, and outlined the procedure which is straightforward to implement in order to obtain 
nongaussianity at an arbitrary order. The results were applied for the case of radiation for illustration, and they are 
easily applied to any other fluid, given nothing but its thermodynamic properties specified by the equation of state, 
or equivalently, the temperature dependence of its energy density. Fundamentally, these follow from the partition 
function. 

Another question is whether there are realistic cosmological models in which statistical thermal fluctuations are 
dominant instead of the usual quantum fluctuations to seed the large scale structures? We believe the cyclic in- 
flationary scenario presents a plausible framework where this indeed turns out to be the case, and in a companion 
paper we shall apply the results obtained here to study this scenario in detail and show that there are parameter 
regions compatible with the present observations and falsifiablc predictions for both the tcnsor-to-scalar ratio and 
nongaussianity. 



TB would like to thank Marco Peloso for raising some important issues about the nature and evolution of thermal 
fluctuations. TB's research has been supported by the LEQSF(2011-13)-RD-A21 grant from the Louisiana Board of 
Regents. TK is supported by the Research Council of Norway, AM is supported by the Lancaster-Manchester-Shefheld 
Consortium for Fundamental Physics under STFC grant ST/J000418/1. 



Here we are going to take a closer look at how the fluctuations in energy in a given physical volume is converted 
to the power spectrum. To begin with we note that the thermodynamic computation of energy fluctuations is done 
in an "adiabatic" approximation scheme where we ignore the cosmological evolution. Thus tells how in a given 
"Eucledian" time slice the "relative energy" in a given physical volume fluctuates. The question we are interested 
in asking is that for a given time slice, if we know how Sl depends on L, how can we compute the Fourier components 
of relative density fluctuations in that same time slice. To understand such a "kinematical/statistical" relationship 
we can therefore work with physical coordinates and momenta. Once we derive the relationship, it is relatively easy 
to convert the results in "comoving" language which is more useful for cosmology. 

To keep things finite we first choose a "fiducial" volume in our universe, V, which is big enough that it contains all 
the relevant scales we are interested in i.e. , L <^ L. We are not going to provide all the rigorous details/justifications 
for doing this, which arc presented in [49j . We are also going to assume periodic boundary conditions, so that physical 
momenta's are confined to integral values: 



V. 



CONCLUSIONS 
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Appendix A: Window Function 




n . 



(Al) 



We can now define the Fourier components via 




(A2) 



(A3) 
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and the relative density fluctuations are defined in the usual way 

d{x) 



- where (p) ^ ^^^^^ '^'^ 



(A4) 



(P> " ' V 

Now, what we are interested is to obtain a statistical measure of S. However, statistically one expects 

{S)=0. (A5) 

Thus, one has to obtain an estimate of (5k using a root mean square approach. The standard approach is to use the 
correlation function: 



^{SnSn' cxp 



i27rn • (x + r) 



C(r)^(5(x)<5(x + r))^ 

^(SnSn' cxp 



(A6) 





i27rn' • x 


exp 






L 



i27r(n — n') • x 





i27rn • r 


exp 






I 



exp 



i27rn • r 



Our final step is to go from the discretized momentum space to a continuum limit {V ^ oo limit). To achieve this, 
as is always done in statistical mechanics, wc have to convert sums into integrals and go from the n-space to k-space: 



d'n{\Snn = / d'k Si 



V 



^"''^ (2^ 

since the n and k spaces arc related via the density of states: 



V 



(l^nl 



■d-^k. 



Thus the correlation function can now be expressed as 

z27rn • r 



L 



(27r)3 



d k(5k exp(-ik • r) , 



(A7) 



(A8) 



(A9) 



In other words, 5^ is just the fourier transform of the correlation function. Since the correlation function is a physical 
quantity i.e. , it does not depend on the fiducial volume over which the averaging is performed, it is also clear that 
(5^ is also a physical quantity. 

Using similar analysis one can also statistically compute the mass / energy variance in a given physical volume. We 
find (the derivation is straight forward and given in section 13.3 of [49!|) 



Sl = J2^'r.W\27rnL/L), 



where the window function is defined as 



W[kL] 



J exp[zk-y] d^y 



(AlO) 



(All) 



Again we can pass from the discrete n-space to the continuum momentum /c-spacc to find 



Si 



(2' 



^ J d^k SlW^ikL) J dk k^'SlW^kL) = ^ J dk k'P{k)W^{kL) (A12) 
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One can write down an analytical expression for the window function (|A11[) in terms of Bessel functions. However, 
to avoid some technical complications in literature often other window functions are used such as 

iyG(kL) = exphi(kL)2]. (A13) 

We emphasize that in the above formulas k and L represent physical momenta and lengths respectively. The inversion 
thus gives us 

The final step step involves going from the physical fourier components to the comoving fourier modes. The two 
are defined via 

1 

(2V)3 



where we now denote the physical coordinates as vectors x = a(t)x. It is clear now that 

<i>l = a-^ml. (A17) 



Appendix B: Quantum vs. thermal fluctuations 

To estimate the relative contributions to the density perturbations originating from thermal statistical vs. the 
quantum vacuum fluctuations, let us consider an ideal fluid with a constant equation of state parameter < w < 1. 
This corresponds to having a polynomial pressure as a function of temperature: 

.4 fT 



p(T) — m — and p ~ wp (Bl) 

\m J 

where m is a mass scale associated with the fluid. Thus according to our formula for the spectrum, it's parametric 
dependence on the temperature is given by 

In contrast, the quantum vacuum fluctuations yield a power-spectrum that is proportional to the density of the 
background fluid [l|: 

We immediately notice that the vacuum fluctuations have an extra Planck suppression: 



(B5) 



For instance, for radiation, since if = 1/3, we have 



More generally, for the range 1 > w > 1/3, since both the exponents in (jB4[) are positive, as long as T, m < Mp, the 
vacuum fluctuations are suppressed as compared to the random thermal fluctuations. If 1/3 > if > 0, then depending 
upon the value of w, if T is sufficiently larger than m, the vacuum fluctuations may be able to dominate over the 
thermal fluctuations. However, in most physical scenarios one expects T < m for the validity of the physics involving 
the thermal fluid. Thus in most physical scenarios we actually expect the thermal fluctuations to dominate the show, 
but the formal condition for this to happen is given by 
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Appendix C: Pressure perturbation 



For a fluid in thermal equilibrium, the thermodynamic properties of the fluid determine the density fluctuation. 
This, through Einstein's field equations, is related to the gravitational potential which in turn determines the pressure 
fluctuation. We can thus compute the explicit form of the latter. This is a straightforward but tedious task and we 
omit the details of the algebra herc^°. The result is 



For radiation this yields = (2/3)(aT^/A:Mp)^. This is certainly different from c^^^, which again highlights the 
difference between thermal and hydrodynamic perturbations. In particular, the properties of the fluid and the back- 
ground expansion need to be taken into account (both of which in principle is determined by the temperature) , but in 
addition the relation is scale dependent. In fact, approaching the small-volume limit /c — )■ oo the pressure fluctuation 
becomes negligible, and at the largest-scale limit fc — >■ the result formally diverges, but of course thermal correlations 
are only expected to exist as long as fc/a < H . 
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